Current Density and Reversible Magnetization of HgBa 2 Ca 2 Cu30 x Superconductors 
Containing Randomly Oriented Columnar Defects 
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Bulk poly crystalline HgBa2Ca2Cu30a; materials were irradiated with 0.8 GeV protons to form ran- 
domly oriented columnar defects, by induced fission of Hg-nuclei. Proton fluences from to 35 x 10 16 
cm -2 were used to install defects with area densities up to a "matching field" of 3.4 Tesla. Studies 
were conducted on the dependence of the equilibrium magnetization and the intragrain persistent 
current density on temperature and applied magnetic field, at various defect densities. The magne- 
tization was modeled using London theory with the addition of vortex-defect interactions. 
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I. INTRODUCTION 

For strong pinning of vortices in high-T c supercon- 
ductors, some form of correlated disorder generally is 
most effective. Most widely studied have been colum- 
nar defects, which typically are formed by irradiation 
with energetic heavy ions; such particles are highly ioniz- 
ing and create tracks of amorphous material along their 
path.EJ Many morphologies have been investigated, rang- 
ing from the familiar parallel tracks to inclined defects, 
crossed arrays with "designer splay," etc. One problem 
with heavy ions is their limited range, typically a few 
tens of micrometers. To overcome the limited range of 
heavy ions, Krusin-Elbaum et al.0 demonstrated an in- 
direct formation of columnar defects using deeply pen- 
etrating 0.8 GeV protons. In this process, an incident 
proton is absorbed by a heavy nucleus (Hg) within the 
material and excites it into a highly energetic state. As a 
consequence, the nucleus breaks into two particles with 
similar mass, each having ~ 100 MeV energy. The fis- 
sion fragments form randomly oriented columnar defects 
(CD's) deep within the superconductor. In a different ap- 
proach to create pins deep inside a material, Weinstein 
and coworkers have created CD's by introducing 235 U as 
a dopant into the superconductor, which is-induced to 
fission by irradiation with thermal neutrons .□ 

Vortex pinning by proton-generated fission-defects was 
first demonstrated in Bi-cuprate materials £lul In fact, 
the proton-based fission process is generally applicable 
in cuprate superconductors that contain a sufficient den- 
sity of heavy nucleifi High-T c materials containing ran- 
domly oriented CD's exhibit a variety of interesting phys- 
ical phenomena, e.g., temperature independent quantum 
tunneling of vortices in Bi-2212 materials. El In studies of a 



series of Hg-cuprates^containing 1, 2, and 3 adjacent CuO 
layers, it was shownQ that sufficiently high superconduc- 
tive anisotropy can lead to a rescaling of the splayed land- 
scape of random CD's. As a consequence, the pinning ar- 
ray and the applied field are 'refocused' toward the crys- 
talline c-axis, even in polycrystalline materials. This is a 
useful feature that we will use when interpreting present 
studies of the equilibrium magnetization. 

In this paper, we describe how the addition of ran- 
dom CD's changes the superconductive properties of 
HgBa2Ca2Cu30 :! ; superconductors. Specific topics in- 
clude the intragrain persistent current density </, as ob- 
tained from the magnetization; establishment of the op- 
timum defect density; and analysis of the changes in the 
mixed state reversible (or equilibrium) magnetization re- 
sulting from interactions between vortices and randomly 
oriented columnar defects. 



II. EXPERIMENTAL ASPECTS 

Samples for study were bulk polycrystalline 
HgBa2Ca2Cu30a; materials (Hg-1223) containing sets of 
3 adjacent oxygen-copper layers. Small pieces, typically 
30 mg mass and ~ 1 millimeter thickness, were all cut 
from the same pellet. The samples were irradiated at 
room temperature in air with 0.8 GeV protons at the Los 
Alamos National Laboratory. Proton fluences <fi p were 0, 
1.0, 3.2, 10, 19, and 35 x 10 16 protons/cm 2 , as determined 
from the activation of Al dosimetry foils. The resulting 
density of fission events is N/V = (ppafUHg- Here nn g 
is the number density of Hg nuclei and <jf 80 mil- 
libarns is the cross-section for inducing a prompt fission 
of Hg nuclei. Each fission produces one CD, so N/V is 
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the volume density of defects, also. One can convert this 
into an approximate area density of CD's by multiply- 
ing by the track length ~ 8 /im, and the area density 
can be reexpressed in units of a matching field i?$ by 
multiplying by the flux quantum <i>o. Resulting values 
for the defect density are B$, = 0, 0.1, 0.3, 1.0, 1.9, and 
3.4 Tesla, respectively. The fission process is random 
in direction, giving randomly oriented CD's. Transmis- 
sion electron microscopy (TEM) studies of both Bi-2212 
and Hg-cupratesB have demonstrated the presence of ran- 
domly oriented columnar defects in these high-T c mate- 
rials. 

The superconductive properties of the virgin and 
processed samples were investigated magnetically. A 
SQUID-based magnetometer (model Quantum Design 
MPMS-7), equipped with a high homogeneity 7 T super- 
conductive magnet, was used for studies in the temper- 
ature range 5-295 K, in applied fields to 6.5 T. The su- 
perconductive transition temperatures T c were measured 
in an applied field of 4 Oe (0.4 mT) in zero-field-cooling 
(ZFC) and field-cooling (FC) modes. The resulting val- 
ues for the onset temperature T c are 133, 134, 132.3, 132, 
131, and 129.3 K, respectively. Values for the Meissner 
(FC) fraction ~4ttM/H lie in the range of 40- 50%, ex- 
cept at the highest fluence where the fractional flux ex- 
pulsion was 29%. 

The isothermal magnetization M was measured as a 
function of applied magnetic field. Below T c and be- 
low the irreversibility line, the magnetization was hys- 
teretic due to the presence of intragrain persistent cur- 
rents. From the magnetic irreversibility AM = [M(H J. 
) — M(H t)], the persistent current density J was ob- 
tained using the Bean critical state relation J cx AM/r, 
where r « 4 /jm is the mean grain radius. Measure- 
ments of the background magnetization in the normal 
state above T c were used to correct the data in the su- 
perconducting state, in order to obtain the equilibrium 
magnetization M eq . Above the irreversibility line where 
AM = 0, this process yields M eq directly; near the 
irreversibility line where AM is small, we obtain M eq 
from the (background-corrected) average magnetization 
[M(H I) + M(H |)]/2, as illustrated later. 

III. EXPERIMENTAL RESULTS 

We begin with a discussion of the effects of randomly 
oriented columnar defects on the irreversible properties of 
the Hg-1223 superconducting materials. The addition of 
correlated disorder increases the pinning of vortices, of- 
ten quite significantly. An example is contained in Fig. 1, 
which shows the magnetization M(H) at 60 K vs magne- 
tizing field H, for various defect densities. With increas- 
ing fluence, the "hysteresis loops" increase in width and 
become symmetric about the M = axis. (The asym- 
metry for the virgin sample, coupled with the fact that 
the decreasing field branch lies near the M = axis, sug- 



gests that surface barrierstHlJ contribute to the observed 
hysteresis in this case. It is interesting to note that a rel- 
atively low dosage of CD's, Bg, =0.1 T, symmetrizes the 
M{H) loop significantly; in particular, the magnitude of 
M in the decreasing field branch is much larger. These 
features imply that the addition of CD's suppresses the 
surface haxrier, as recently discussed by Koshelev and 
Vinokur£3) 

Using the Bean model, one may obtain the persistent 
current density J(H,T). For these weakly linked poly- 
crystalline materials, the magnetization reflects the in- 
tragranular current density. Some results of this analysis 
are shown in Fig. 2 and Fig. 3 as plots of J versus 
temperature T, in applied fields of 0.1 and 1 Tesla, re- 
spectively. The enhancement in J is modest in low fields; 
this can be attributed to the presence in the virgin sam- 
ple of both naturally occurring defects that provide some 
pinning and the likely influence of surface barrier effects, 
as already noted. In higher fields, the contribution of 
the random CD's is more apparent and J is enhanced by 
about an order-of-magnitude. For high-T c materials with 
these angularly dispersed defects, one generally achieves 
the maximum J at some defect density _B$ near 0.2-2 T. 
This optimization is shown by the insets in Figs. 2 and 
3. For an applied field of 0.1 T (Fig. 2 inset), J at 
100 K is largest for B$ = 0.1-0.3 T. For the second ex- 
ample with a field of 1 T (Fig. 3 inset), J at 60 K is 
largest for £?$ near 2 T. Qualitatively, the optimum de- 
fect density in each case is comparable with the vortex 
density, i.e., about 2x larger. At higher defect densities, 
J decreases. This may be due in part to the presence of 
additional CD's that help to initiate the hopping of vor- 
tices to nearby empty pinning sites. A second, potential 
influence is a suppression of T c ; however, the change in 
T c is small and these measurements are far from T c , so 
this contribution is expected to be small. A third and 
sizable effect is an irradiation-induced suppression of the 
vortex line energy and pinning energy, due to changes in 
the London penetration depth; this will be discussed in a 
following section. The optimum defect density depends 
weakly on the temperature as well. 

Next we consider the equilibrium magnetization in 
the mixed state. Figure 4 illustrates the method used 
to obtain M eq by plotting the background-corrected ex- 
perimental magnetization for the virgin material (open 
squares) versus applied magnetic field /j,qH on a loga- 
rithmic scale. The closed symbols show the average M, 
which provides a very good approximation to M eq when 
the hysteresis is small. The solid line is a fit to conven- 
tional London theory, which provides that M eq is directly 
proportional to (1/A^) xln(0 B&fH), where X a b is the in- 
plane London penetration depth, (3 is a constant of order 
unity, and B c i is the upper critical field. t3 The straight 
line fit in Fig. 4 shows that the London ln(-ff) logarithmic 
variation describes the field dependence well. Results for 
M eq in the virgin Hg-1223 material at other temperatures 
are presented in Fig. 5. One sees from the linear regres- 
sion lines that the London field dependence is followed 
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over a wide range of fields and temperatures. Deviations 
from linearity occur at low fields when H approaches H c \ 
where simple London theory is not valid, and at low tem- 
peratures where the materials are sufficiently hysteretic 
that the averaging procedure illustrated in Fig. 4 is not 
valid. The slopes of the curves in Fig. 5 provide values 
for the London penetration depth \ a b{T) in the virgin 
material. These results will be compared and contrasted 
with values deduced for the irradiated Hg-1223. 

We now address the question of how the addition of 
randomly oriented columnar defects modifies the equilib- 
rium magnetization in the mixed state. The qualitative 
effects of the CD's are illustrated in Fig. 4, which in- 
cludes data for M eq at 77 K for irradiated samples with 
defect densities B<$> rj 1.0 and 3.4 T. Comparing these 
data with the virgin curve shows that the CD's reduce 
the magnitude of the equilibrium magnetization consider- 
ably, especially at lower fields. In addition, they generate 
a pronounced deviation from the "standard" London field 
dependence, giving M eq an "S"-like dependence on field. 
This "anomalous" behavior has been observed previously 
in cuprates containing parallel columnar defects formed 
by 5.8 GeV Pb-kms, in thallium-based single crystal^ 4 !; 
in Bi-2223 tapesB ; and in Bi-2212 single crystals.BB 

These changes in the equilibrium magnetization have 
been attributed to magnetic interactions between the vor- 
tex lattice and the columnar defects. By occupying a 
pinning site, a vortex gains pinning energy. This reduc- 
tion in system energy must exceed the energy increase 
arising from direct intervortex repulsion when a vortex 
is displaced from its natural position in the lattice to a 
particular columnar defect. For a defect geometry with 
parallel tracks thaLhave a Poisson distribution of separa- 
tions, Wahl et al.Efl obtained an expression for M eq that 
describes reasonably well the S-shaped field dependence 
in Tl-cuprate crystals containing parallel CD's. 

In the Hg-cuprates investigated here, one might ex- 
pect that the randomly oriented columnar microstruc- 
ture should entangle the vortices. Consequently it is 
somewhat curious that the M eq in polycrystalline ma- 
terials with random CD's can resemble single crystals 
with parallel defects. We have suggestedta that this 
similarity originates from an anisotropy-induced "refo- 
cussing" of the defects and field toward the crystalline 
c-axis. For highly anisotropic single crystals, it is long 
recognized, that only the normal component of field is 
effective.lla More recent studies of polycrystalline Hg- 
cuprates demonstrated a recovery of vortex variable 
range hopping (VRH), very similar to that observed, 
in YBaCuO single crystals containing parallel CD's.EHl 
Among the Hg-cuprates with 1, 2, and 3 adjacent Cu- 
O layers, the recovery of VRH was most pronounced in 
the Hg-1223 material with the largest mass anisotropy 
parameter 7. According to the theoretical development, 
the complexity of randomly oriented columnar defects in 
a polycrystalline material is reduced to some degree by 
a large superconducting anisotropy, restoring the simple 
physics of a crystal with parallel pins. Thus vortex pin- 



ning by random CD's provides a reasonable qualitative 
explanation for the reduction in equilibrium magnetiza- 
tion, despite the complexity of the materials in real space. 

Next we,-iise the vortex-defect interaction theory of 
Wahl ct alM to model our experimental data. From the 
theory, one has 
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where Eq — [<I>o/47rA a {,] 2 is the line energy, J7 is the pin- 
ning energy, and B = (H + AttM) = H since M is small. 
The first term is the conventional London expression. 
The second is an added term to account for interactions; 
it is most significant in intermediate fields and it vanishes 
in large fields B In modeling of data at various 

temperatures T (77, 90, 100, and 110 K), our objective 
was to maintain an internally consistent set of parame- 
ters. The results are shown as solid lines in Fig. 6a for 
the sample with B<j, ~ 1 T and in Fig. 6b for the sample 
with a higher defect density. Given the complexity of 
the polycrystalline material and of the vortex and defect 
arrays, the description of the experimental data (filled 
symbols) is relatively good. The values for the pinning 
energy are a significant fraction of the line energy and 
are quite reasonable, with Uq = 0.65eo an d 0.8eo ; re- 
spectively. To obtain reasonable modeling of the data 
for the two irradiated materials, we used values of 2.0 
and 2.5 T, respectively, for the effective defect densities 
Bg, . The differences from the nominal values (calculated 
from the proton fluences) may arise from some overlap of 
tracks at the highest fluence combined with uncertainty 
in calculating the production rate for CD's. It is also 
possible that the effective values for i?$ compensate for 
other factors, such as residual entanglement and angular 
distribution in the CD array. Indeed, for more highly 
anisotropic Tl-2212 materials (7 > lOO—compared with 
7 ~ 60 in Hg-1223), similar modelings was achieved 
with B$ values within 20% of those calculated from the 
proton fluence. 

An important feature of any superconductor is the 
London penetration depth. Figure 7 summarizes the re- 
sults for these materials in a plot of 1/A 2 h versus tem- 
perature. Values for the virgin sample were obtained 
from a standard London analysis of the data in Fig. 5. 
An extrapolation of the Ginzburg-Landau linear depen- 
dence near T c (dashed line) to T = yields a value 
A(T = 0) = 174 pp, which is comparable with ear- 
lier determinations l3 The values for the materials with 
CD's were obtained from the modeling procedure. All of 
these results are reasonably behaved, with 1/A 2 b vary- 
ing linearly with T at high temperature. It is clear 
that proton irradiation increased A significantly. This 
is expected, from two theoretical perspectives. First, 
conventional GLAG theoryt 2 ] predicts that the penetra- 
tion depth A 2 increases as (1 +£0/^) when the electronic 
mean free path i is reduced by scattering, as produced 
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by irradiation-generated defects from neutrons and sec- 
ondary protons released by spallation. Second, the the- 
ory of Wahl-Buzdinli-3 provides that introducing CD's in- 
creases the penetration depth as 

A" 2 (B$) = A" 2 (S$ = 0) x [1 - 27ri? 2 S<i,/$o] (2) 

where R is the radius of the columnar track. The com- 
bination of these two effects leads to the significant, pro- 
gressive increases in A observed in Fig. 7. To compare 
with Eq. 2, the inset of Fig. 7 shows 1/A 2 , linearly ex- 
trapolated to T = 0, plotted versus defect density. The 
open symbols denote values for Bq, calculated from the 
proton fluence, while the closed symbols are values used 
for the modeling, as discussed above. The dashed line 
in the inset illustrates the theoretical dependence in Eq. 
2 and shows a qualitative agreement with the data; its 
slope corresponds to a value R « 8.4 nm, somewhat be- 
low the 11 nm value,-obtained for Tl-2212 materials. As 
discussed previously)!! it is likely that the CD's have a 
large effective radius due to oblique passage of ion tracks 
through CuO planes and a "halo" ©£, oxygen depletion 
around the central amorphized zone.Eia 

An interesting consequence of a larger London pene- 
tration depth is that the vortex line energy e decreases 
significantly. As a result, one can expect the vortex pin- 
ning energy of a CD and its pinning force to decrease, re- 
ducing its effectiveness. This effect counteracts the nom- 
inally positive influence of increasing the defect density 
£>$ . It is likely that these competing effects play a major 
role in determining the range of -B$ that maximizes the 
current density J, as illustrated in the insets of Figs. 2 
and 3. This competition is reminiscent of the formation 
of (point-like) defects in YBa2Cu307_,5, where progres- 
sive removal of oxygen creates defects, but also reduces 
their effectiveness, through increases in the superconduc- 
tive length scales£ 



IV. SUMMARY 

We irradiated polycrystalline HgBa2Ca2Cu30a; mate- 
rials with 0.8 GeV protons to produce randomly oriented 
columnar defects via a fission process. Enhancements in 
vortex pinning increased the persistent current density, 
with the optimum defect density depending on the range 
of field and temperature. Adding random linear defects 
significantly reduced the magnitude of the equilibrium 
magnetization and changed its dependence on magnetic 
field from a simple London ln(iJ) form (as observed for 
the virgin material) to a more complex "S-shaped" de- 
pendence. Invoking an anisotropy-induced "refocusing" 
of the vortex-defect array, we have modeled these re- 
sults using the theory of Wahl-Buzdin that incorporates 
vortex-defect interactions. This analysis shows that the 
addition of random columnar defects increases the Lon- 
don penetration depth markedly, which diminishes their 
effectiveness for pinning vortices. Overall, however, the 



formation of randomly oriented columnar defects both 
enhances the current-carrying performance of the mate- 
rial and creates interesting superconductive systems for 
study of the interaction of vortices with correlated disor- 
der. 
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FIGURE CAPTIONS 

Fig. 1. The magnetization of polycrystalline Hg-1223 
materials at 60 K, versus applied magnetic field. Sam- 
ples were irradiated with various flucnces <j) p of 0.8 GeV 
protons, as shown, to create randomly oriented columnar 
defects with approximate "matching fields" 

Fig. 2. The intra-granular persistent current density 
J vs temperature for various irradiated materials, mea- 
sured in applied magnetic field of 0.1 T. The inset shows 
that J at 100 K is maximized at a defect density of 0.1- 
0.3 T. 

Fig. 3. The intra-granular persistent current density 
J vs temperature for various irradiated materials, mea- 
sured in applied magnetic field of 1 T. Inset: J at 60 K 
versus defect density £?$ exhibits a maximum at a defect 
density of -2 T. 

Fig. 4. The magnetization of virgin Hg-1223 at 77 K, 
versus applied magnetic field (log scale). Open squares 
show the measured M in increasing and decreasing field 
history, while closed squares show the average M k, M eq . 
The solid line is a fit to the conventional London \n(H) 
dependence. Data for Hg-1223 containing randomly ori- 
ented columnar defects are included for comparison. 

Fig. 5. The equilibrium magnetization M eq (H,T) 
plotted vs ln(H) for unirradiated Hg-1223 at the tem- 



peratures shown. Lines are fits to conventional London 
theory. 

Fig. 6. Plots of M eq (H,T) vs. \n(H) for irradiated 
Hg-1223 with (a) B$ ~ 1 T and (b) with ~ 3.4 T. 
Symbols are experimental results, while lines show mod- 
eling of data using Eq. 1 where effects of vortex-defect 
interactions are included. 

Fig. 7. The temperature dependence of the London 
penetration depth 1/A 2 . Values for the virgin sample 
come from a conventional London analysis of the data 
in Fig. 5; for the irradiated materials, values come from 
modeling the equilibrium magnetization shown in Fig. 6. 
Inset shows values for 1/A 2 , extrapolated to T = 0, plot- 
ted versus defect density; open symbols denote -B$ values 
calculated from the proton fluence and closed symbols are 
values used in the modeling. 
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